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Abstract — We present a family of non-additive quantum 
codes based on Goethals and Preparata codes witli parameters 
((2™, 22"-5'"+\ 8)). Tlie dimension of these codes is eiglit times 
higher than the dimension of the best known additive quantum 
codes of equal length and minimum distance. 

Index Terms — Non-additive quantum code, Goethals code, 
Preparata code 

I. Introduction 

Most of the known quantum error-correcting codes 
(QECCs) are based on the so-called stabilizer formalism which 
relates quantum codes to certain additive codes over GF{4) 
(see, e. g., f3l, f7|). It is known that non-additive QECCs can 
have a higher dimension compared to additive QECCs with the 
same length and minimum distance ||5), |[T4|, |[T7|, |[T8|. All 
these examples of non-additive QECCs are examples of so- 
called codeword stabilized quantum codes which are obtained 
as the complex span of some so-called stabilizer states, which 
correspond to self-dual additive codes. In |9| we have extended 
the framework of stabilizer codes to the union of stabilizer 
codes (see jS)). This allows to construct non-additive codes 
from any stabilizer code. In general, these non-additive QECCs 
correspond to non-additive codes over GF{4) which can be 
decomposed into cosets of an additive code which contains 
its dual. Using a construction similar to that of so-called CSS 
codes (see Q, |[T5|), families of non-additive quantum codes 
based on the binary Goethals and Preparata codes were derived 
in Q. Here we present a new family of non-additives quantum 
codes which have a dimension that is eight times higher than 
the dimension of the best known additive quantum codes. 

II. Union Stabilizer Codes 

A. Stabilizer codes 

We start with a brief review of the stabilizer formalism for 
quantum error-correcting codes and the connection to additive 
codes over GF{A) (see, e.g., ||3], ||7|). A stabilizer code 
encoding k qubits into n qubits having minimum distance d, 
denoted by C = [[ri, fc,c?]], is a subspace of dimension 2''" of 
the complex Hilbert space (C^)**" of dimension 2". The code 
is the joint eigenspace of a set of n — k commuting operators 
Si, ... , Sn~k which are tensor products of the Pauli matrices 



or identity. The operators Si generate an Abelian group S 
with 2"^''' elements, called the stabilizer of the code. It is 
a subgroup of the n-qubit Pauli group Vn which itself is 
generated by the tensor product of n Pauli matrices and 
identity. We further require that S does not contain any non- 
trivial multiple of identity. The normalizer of S in Vn, denoted 
by TV, acts on the code C = [[n, fc, d]]. It is possible to identify 
2fc logical operators Xi, . . . , Xk and Zi, . . . ,Zk such that 
these operators commute with any element in the stabilizer S, 
and such that together with S they generate the normaUzer 
Af of the code. The operators Xi mutually commute, and so 
do the operators Zj. The operator Xi anti-commutes with the 
operator Zj if i — j and otherwise commutes with it. 

It has been shown that the 7i-qubit Pauli group corresponds 
to a symplectic geometry, and that one can reduce the problem 
of constructing stabilizer codes to finding additive codes over 
GF{4) that are self-orthogonal with respect to a symplectic 
inner product Q, Up to a scalar multiple, the elements 
of Vi can be expressed as a^cr^ where (a, b) e is a 
binary vector. Choosing the basis of Gi^(4), where 



w is a primitive element of GF(4) with u 



1 = 0, we 



get the following correspondence between the Pauli matrices, 
elements of GF(4), and binary vectors of length two: 



operator 


GF(4) 
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(10) 
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(01) 



This mapping extends naturally to tensor products of n Pauli 
matrices being mapped to vectors of length n over GF(4) or 
binary vectors of length 2n. We rearrange the latter in such a 
way that the first n coordinates correspond to the exponents 
of the operators ax and write the vector as {a\b), i.e.. 
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{a\h) = {g^\g^). (1) 
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Two operators corresponding to the binary vectors (a|6) and 
(cjd) commute if and only if the symplectic inner product 
a ■ d — b ■ c = 0. In terms of the binary representation, 
the stabilizer corresponds to a binary code G which is self- 
orthogonal with respect to this symplectic inner product, and 
the normalizer corresponds to the symplectic dual code G*. 
In terms of the correspondence to vectors over GF(4), the 
stabilizer and normalizer correspond to an additive code over 



GF{4) and its dual with respect to an symplectic inner 
product, respectively, which we will also denote by C and C* . 
The term additive quantum code refers to this correspondence. 
The minimum distance d of the quantum code is given as the 
minimum weight in the set C* \ C C GF{A)'^ which is lower 
bounded by the minimum distance d* of the additive code C* . 
If d — d* , the code is said to be pure, and for d > d* , the 
code is said to be pure up to d* . 

Fixing the logical operators Xi and Zj, there is a canonical 
basis for the additive quantum code C. The stabilizer group 
S of the quantum code together with the logical operators 
Zj generate an Abelian group of order 2" which corresponds 
to a self-dual additive code. The joint +l-eigenspace is one- 
dimensional, hence there is a unique quantum state |00 . . . 0) G 
C stabilized by all elements of S. An orthonormal basis of the 
code C is given by the states 

\ii^2 ■■■ik) = ---^k lOOTTTO) , (2) 
where (11^2 • ■ ■ ifc) G IFl- 

B. Union stabilizer codes 

The stabilizer group S gives rise to an orthogonal decompo- 
sition of the space (C^)*^" into common eigenspaces of equal 
dimension. The stabilizer code C is the joint +l-eigenspace 
of dimension 2*^. In general, the joint eigenspaces of S can 
be labeled by the eigenvalues of a set of n — fc generators of 
S. Moreover, the n-qubit Pauli group Vn operates transitively 
on the eigenspaces. Hence one can identify a set T C Vn of 
2n~k operators such that 

(C2)®" = 0tC. (3) 

tGT 

Note that each of the spaces tC is a quantum error-correcting 
code with the same parameters as the code C and stabilizer 
group tSt^^. The decomposition (jsjl corresponds to the de- 
composition of the n-qubit Pauli group Vn into cosets with 
respect to the normalizer M of the code C and likewise to the 
decomposition of the full vector space G-F(4)" into cosets of 
the additive code C* . 

The main idea of union stabilizer codes is to find a subset 
To of the translations T such that the space ©fg^-^ tC is a 
good quantum code (see ||8|, |^). 

Definition 1 (union stabilizer code): Let Co = [['^j^]] a 
stabilizer code and let To — {ti, . . . ,^^1 be a subset of the 
coset representatives of the normalizer Aq of the code Co in 
Vn- Then the union stabilizer code is defined as 

C - tCo. 

teTo 

Without loss of generality we assume that To contains identity. 
The dimension of C is K2^, and we will use the notation 

C = {{n,K2\d)). 

Similar to ([2]) a canonical basis of the union stabilizer code C 
is given by 

\3\m2-.-ik) = tjX\' ---Xl' lOOTTO) , (4) 



where j = 1, . . . ,K, {iii2 ■ ■ - ik) G 1^2' logical 
operators of the stabilizer code Co. 

In order to compute the minimum distance of this code, we 
first consider the distance between two spaces tiCo and t2Co- 
As for a fixed stabilizer code Co two spaces tiCo and i2Co are 
either identical or orthogonal, we can define the distance of 
them as follows: 

dist(iiCo,t2Co) := min{wgt(p): p^Vn\ phCa = <2Co}. 

(5) 

Here wgt(p) is the number of tensor factors in the n-qubit 
Pauli operator p that are different from identity. Clearly, 
dist(tiCo, i2Co) = dist(t^^tiCo,Co). The two spaces are 
identical if and only if t^^ti is an element of the normalizer 
group Ao, or equivalently, if the cosets Co +ti and C^+t2 of 
the additive normalizer code Co are identical. (Note that we 
denote both an n-qubit Pauli operator and the corresponding 
vector over GF(4) by ti.) Hence the distance (jsjl can also be 
expressed in terms of the associated vectors over GF(4). 

Lemma 2: The distance of the spaces iiCo and t2Co equals 
the minimum weight in the coset Cq + ti — ^2- 
Proof: Direct computation shows 

dist(iiCo, tsCo) = dist(Co* + ti, Cq* + ^2) 
= dist(Co* + (ti -t2),Co*) 
= min{wgt(c + - ^2) : c € Cq} 
= min{wgt(u) : u G Co + ti - 12}. ^ 

While the distance between the cosets Cq + tj is an upper 
bound on the minimum distance of the union code C, the true 
minimum distance can be derived from the following code 
over CF(4). 

Definition 3 (union normalizer code): With the union sta- 
bilizer code C we associate the (in general non-additive) union 
normalizer code given by 

C* = \J C*+t = {c + tf. ce Co*, j^l,...,K}, 

teTo 

where Cq denotes the additive code associated with the 
normalizer Aq of the stabilizer code Co. We will refer to 
both, the vectors ti and the corresponding unitary operators, 
as translations. 

Theorem 4: The minimum distance of a union stabilizer 
code with union normalizer code C* is given by 

d ^ min{wgt(w) : u e (C* - C*) \ Co} 
> d„,i„(C*) 

= min{dist(c + ti, c' + ti') : ti, ti' S %, c, c' £ Cq 

C+t,:^ C' + ti,}, 

where C* — C* := {a — b: a,b E C*} denotes the set of all 
differences of vectors in C*, and Co < Co is the additive code 
that corresponds to all elements of the stabilizer group S that 
commute with all tj G Tq. 

Proof: Let E E Vn be an n-qubit Pauli error of weight 
< wgt{E) < d. For two canonical basis states jV-'a) and 



as given in Q we consider the inner product 



Definition 6 (Goethals code |7]?j; The Goethals code Q{ra) 
of length 2™ consists of the codewords described by all pairs 
(X, Y) satisfying: 



-^{m...Q\Xl ■■■X^t,Etj,X^ ■■■X^\m...Q) a) \X\ is even, |r| is even, 



If -E e 5 commutes with all tj E Tq, then {tpa\E\ipii) = Sab- 
Otherwise, E ^ C* — C* since < wgt(£') < d, and hence 
the inner product vanishes. ■ 

III. The Binary Goethals and Preparata Codes 

In this section we recall some properties of the binary 
Goethals codes |6) and the Preparata codes [13|. It has been 
shown that variations of these codes have a simple description 



as Z4-linear codes |10|, but in our context the description in 
terms of cosets of linear binary codes is used. 

In the following m is an even integer (m > 4) and 
n — 2™~^ — 1. Leta be a primitive element of the finite field 
GF(2™^^). By iJ,i{z) we denote the minimal polynomial of 
a* over GF{2), i. e., the polynomial with roots for j — i2^. 
The idempotent 9i{z) is the unique polynomial satisfying 

Oiia') = 1 and 9,{a^) = for j ^ 12''. 

Codewords of a cyclic code can be represented by polynomials 
f{z), and we use {f{z); /(I)) to denote the codeword of the 
extended cyclic code obtained by adding an overall parity 
check. Similar, we use {f{z);f{l);g{z);g{l)) to denote the 
juxtaposition of codewords of two extended cyclic codes. 

Definition 5 (Goethals code f^j: The Goethals code G{m) 
of length 2™ is the union of 2™"^ cosets of the linear 
binary code Cg = [2™, 2" - 4m + 2,8]. The code Cg is 
obtained via the |it|ii + w| construction applied to the extended 
cyclic codes Ci and C2. The cyclic code Ci is a single-error 
correcting code with generator polynomial fii{z), and C2 is 
generated by iii{z)iir{z)iJ,s{z) where r — 1 + 2™/^"^ and 
s = 1 + 2™/^^!. The non-zero coset representatives are given 
by {z';l;z'9i{z);0) for i = 1, . . . , n - 1. 

An alternative description of Goethals codes has been given 
in ([1]. The codewords are described by pairs {X, Y) of subsets 
of GF{2™^^). The corresponding codeword is given by the 
juxtaposition of the characteristic functions xx and xy of the 
two set X and Y, i. e. 

{X,Y) = (lx(a');lx(0);ly(a^);ly(0)), 

where is a short-hand for the vector 

Ix(a^) = Ix(a'), . . . , UK"')) 

and 



1 ifxeS, 

if x(is. 



The non-zero elements of X and Y give rise to the polyno- 
mials fx{z) and fviz) given by 



fsiz) 



E 

i=0 



(6) 



xex 

E 



E^- 

yeY 



xex 



E^ 

\xex 



2^y 

yeY 



d) 



xex 



X 



E 



^1 =Ey'- 

vxex I yeY 
In order to relate the two definitions, we distinguish three 
cases. 

1) X — Y\ Conditions c) and d) imply that 'Yl,xex 2^ — 0. 
This is true for all codewords of the cycUc code gener- 
ated by Adding an overall parity check implies 
Condition a). 

T) X = %: The left hand side of Conditions b), c), 
and d) vanishes, so the solutions for Y correspond 
to an extended cyclic code with generator polynomial 

[Lx{z)[Lr{z)^^{z). 

3) X = {0,a; = a*}: From ^ it follows that /y(a) = 
'YliyeY y- Condition b) holds for the set Y corre- 
sponding to fyiz) = z'^9i{z). The left hand side of 
Conditions c) and d) vanishes, so the solutions for Y 
are elements of the extended cyclic code with generator 
polynomial ^r{z)lJ's{z)- As neither r nor s is a power of 
two, the polynomial 9i{z) and hence fyiz) — z^9i{z) 
vanishes for and a'*, i.e.. Conditions c) and d) hold. 
Finally, all codewords of the Goethals code as given in 
Definition [5] are the juxtaposition of two binary vectors of 
even weight, i.e.. Condition a) holds. Hence any codeword 
given by Definition [5] fulfills the conditions of Definition [6] 
The equivalence of the definitions follows from the fact that 
the codes have equal size. 

Next we consider the definition of Preparata codes similar 
to Definition |6] given in f\\. 

Definition 7 (Preparata code [1]): The extended Preparata 
code V{m) of length 2™ and parameter a consists of the 
codewords described by all pairs (X, Y) satisfying: 

a) 1^1 is even, |F| is even. 



b) 



E 

xex 



yeY 



xex \xex / yeY 

Here cr is a power of two and gcd((T ± 1, n) = 1. 
For <T = 2™/-^^^ and n = 2™^^ — 1 we compute 



-1- 2 



)m/2-l 



± 1 



2/2 



T2 =1, 



showing that gcd((T ± l,n) = 1. Hence for this particular 
choice of cr, the Preparata code of Definition [7] contains 
the Goethals code. What is even more, we can describe the 
Preparata code similar to Definition |5] as the union of cosets 



of the linear binary code C-p which contains the linear binary 
code Cg. 

Definition 8: The extended Preparata code V{ra) of length 
2™ is the union of 2™"^ cosets of the linear binary code 
Cv = [2™,2™ - 3to + 1,6]. The code Cv is obtained via 
the + construction applied to the extended cyclic codes 
Ci and C3. The cyclic code Ci is a single-error correcting 
code with generator polynomial iii{z), and C3 is generated 
by where s = 1 + 2™/^~i. The non-zero coset 

representatives are given by (z'; 1; z^9i{z); 0). 
Comparing Definitions |5] and [8] we see that we can use the very 
same coset representatives to construct the Goethals and the 
Preparata code as union of cosets of the linear binary codes Cg 
and C-p, respectively. Moreover, all codes lie between codes 
that are equivalent to the Reed-Muller codes RM{m — 3, m) 
and RM{m - 2, m) = [2™, 2" - rn - 1, 4] (see fu)). This is 
illustrated by the following diagram: 



[2™, 2"* - TO - 1,4] =RM(to-2,to) 



normalizer of the code C is generated by 



[2™, 2™ - 3to+ 1,6] = Cv 



[2™^2'"-4to + 2,8] = Cg 



r{m) = [j^Cv + t, 



- U, Cg + U 



RM(to - 3, to) 

The components of the codes are summarized as follows: 

Ci: cyclic code generated by /ii(z) 

cyclic code generated by /ii(z)/is(z) 
cyclic code generated by ni{z)iJ,r{z)fJ-s{z) 
r = 1 + 2™/2-2^ s = 1 + 2™/2-i 

+ construction applied to the extended cyclic 
codes Ci and C2 

\u\u + v\ construction applied to the extended cyclic 
codes Ci and C3 
n + 1 coset representatives with 



C3: 
C2: 

Cg: 



C 



ti 



(z^ 1; 2*6*1 (z);0) for i = 0, 
for i = n. 



1, 



J0,...,0) 

IV. The Quantum Goethals-Preparata Codes 

Before presenting the new family of non-additive quantum 
codes, we recall Steane's construction to enlarge the dimension 
of CSS codes. 

Theorem 9 (see f^Wj): Let C = [n, k, d] and C" = [n, k' > 
k + l,d'] be linear binary codes with C"*- < C < C. Then 
there exists an additive quantum code C = [[n, k + k' — n,> 
mm{d,3d' /2)]]. Given a generator matrix G of the code C 
and a generator matrix D of the complement of C in C, the 




where A is a fixed-point free linear transformation. 
As the code Cg contains a code that is isomorphic to the 
Reed-Muller code RM{m - 3, to) it follows that Cg <Cg. 
Hence we can apply Steane's construction | [T6[ to the chain 
Cg < Cg < C-p of linear binary codes and obtain an additive 
quantum code with parameters Cq — [[2'" , 2™ — 7m + 3,8]]. 
In a second step we use the K = 2™^^ coset represen- 
tatives ti of the decomposition of both the Goethals and 
the Preparata code. This yields a non-additive code with 
dimension k'^2'^'^-'^"'+^ = 2^ where ^ = 2™ - 5to + 1. 
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Fig. 1 . Structure of the non-additive union normalizer code of the quantum 
Goethals-Preparata codes. 



Theorem 10: Let Co = [[2™, 2"-7to+3, 8]] be the additive 
quantum code obtained from the chain of linear binary codes 
Cg < Cg < C-p using Steane's enlargement construction. 
Furthermore, let Tq = {it,\t^): i,j = 0, 2"^i - 1} where 
ti are the coset representatives used to obtain the Goethals and 
Preparata code. Then the quantum Goethals-Preparata code 
is a union stabilizer code given by Cq and Tq. The minimum 
distance of the quantum Goethals-Preparata code is eight. 

Proof: Let G denote a generator matrix of the code Cg 
and let D be such that ( ^ ) generates C-p. The structure of the 
non-additive union-normalizer code of the quantum Goethals- 
Preparata codes is illustrated in Fig. [T] A generator matrix of 
the normalizer of the additive quantum code Cq is given above 
the horizontal line, while the set of translations is listed below 
the horizontal line. Every codeword of the non-additive union 
normalizer code is of the form 

9 = (5^1.9^) = (ci + V + t,\c2 + W + tj), 

where ci,C2 e Cg ^ [2™, 2™-4m+2, 8] andt;,^ e Cv/Cg. 
For g, g' G C* ,g g' we compute 

dist(g,.g') =dist((ci + w + ti\c2 + w + tj), 
(c'i+w' + i^|4 + w' + t;-)) 
= wgt((c'/ + v" + ti- t[\c2 + w" + tj - t'j)), 



where c" — ci — c'l and c'( — ci — c[ are codewords of Cg, 
and v" = V — v', w" — w ~ w' are codewords of C-p /Cg . In 
general, the weight of g = {g'^\g^) is given by 



1 



wgt((5-^|5^))= -(wgt(5-^) + wgt(g^) + wgt(<?^ +g^)). 
Hence we get 



dist(5,g') = 
+ Jwgt(4' + w" + tj -t') 



(7a) 
(7b) 



2 

+ \ wgt (c" + 4 + v" + w" + U~t'^ + ~t)). (7c) 

By Steane's construction the vectors v" and w" are either 
both zero, or both are non-zero and they are different. For 
v" = w" = 0, we can assume without loss of generality that 



the vectors in (7ai and (7bi are both non-zero. The weight 
of these vectors equals the distance between two codewords 
of the Goethals code, so it is at least 8. For v" ^ Q ^ w" 
the terms (7a i and (jTbji equal the distance of two codewords 
of the Preparata code, so they are lower bounded by 6. We 



will show that for v" ^ w", the vector in (7c i is a non-zero 
codeword of the linear code isomorphic to the Reed-Muller 
code RM{m — 2,m), hence its weight is at least 4. For this, 
consider the vectors a = (ai; 02) — c'l + cj,' + v" + w" ^ 
and b = (61; 62) = ti — t'i + tj — t'j- The coset representatives 
are of the form = (z*; 1; z*0i(z); 0), so the second half 
62 of 6 is a codeword of the extended cyclic code generated 
by 9i{z), while 02 is a codeword of the extended cyclic code 
generated by /ii(2). The intersection of the two codes is trivial, 
so a2 — ^2 only if 02 = 62 = 0. Then wgt(6) < 4 while 
wgt(a) > 6 since 7^ a G C-p. Hence a 7^ 6. ■ 
To our best knowledge, the best additive quantum code with 
the same length and minimum distance has dimension 2'" — 
5to — 2. Codes with these parameters can, e. g., be obtained 
by applying Steane's construction to extended primitive BCH 
codes [2'", 2" -2m -1,8] and [2'", 2" -3m -1,6] (see |^16J). 
In the following table we give the parameters of the first codes 
in these families. Additionally, we give the parameters of the 
non-additive quantum codes derived from Goethals codes in 



Goethals 


enlarged BCH 


Goethals-Preparata 


((64,230,8)) 


[[64, 32, 8]] 


((64,235,8)) 


((256,2210,8)) 


[[256,214,8]] 


((256,2217,8)) 


((1024,2906,8)) 


[[1024,972,8]] 


((1024,2975,8)) 



V. Conclusions 

We have constructed some new non-additive quantum codes 
from nested non-linear binary codes which can be decomposed 
into cosets of linear codes which contain their dual. It is 
interesting to find more good non-linear binary or quaternary 
codes with this property. 



Recently, Ling and Sole have constructed some non-additive 
quantum codes from Z4-linear codes using a CSS-like con- 
struction y_2J. So far it is not clear whether the non-additive 
codes presented here can also be put into the framework of 
Z4-linear codes. 
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